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Abstract 

We study low energy shape oscillations of negative parity in the first and second (isomeric) 
minima in actinides. As a main tool we use the phenomenological Woods-Saxon potential with 
a variety of shape deformations. This allows to include a mixing of various multipolarities when 
considering oscillations with a fixed K quantum number. The phonon energies are determined 
either from the collective Hamiltonian with the microscopic-macrocopic energy and cranking mass 
parameters, or from its simplified version with the constant mass parameters. The results for K w = 
0~,1 ~ in the first minima are in a reasonable agreement with experimental data, including predicted 
El transitions; the K n = 2~ energies are systematically overestimated. In the second minimum, 
as compared to the data for 240 Pu and 236 U, our calculated K =1,2 energies are overestimated 
while the K = energies are three or more times too large. This signals either a non-collective 
character of the experimentally assigned K = states or a serious flaw of the model in the second 
minimum. More data on the K = 0, I 11 = 1~ collective states in the second minima of other nuclei 
are necessary to resolve this issue. 
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I. INTRODUCTION 



Recently, a considerable amount of data has been gathered on the nuclear states in the 
second well in the actinide region, especially in 240 Pu, see [l| and references cited therein. 
Various rotational bands, most of them of negative parity, have been identified in this nucleus 
with moments of inertia characteristic of superdeformation (SD). Such spectroscopic data 
provide a much wanted test for nuclear models that were originally fitted in a region of 
deformations around the ground state (g.s.) minima: How much of their predictive power 
is left beyond that region? Negative-parity shape oscillations are natural candidates for 
the observed low energy band- heads at the SD. Here we study them within the much used 
microscopic-macroscopic model based on the deformed Woods-Saxon potential. The found 
properties of excitations in the shape-isomeric minima provide both a prediction and a test of 
the model. We check the used method by first applying it to the negative-parity excitations 
in the first minima in actinides, on which data are more abundant. 

Discovery of fast fissioning states in actinides by Polikanov et al. 2] was interpreted 
soon afterwards by Flerov and Druin [3] and Strutinsky [4] as evidence for very deformed 
secondary minima in these nuclei. Later experiments provided support for this guess, with 
measurements of rotational band in 240 Pu [5] and then, of the quadrupole moment in the 
second well in 239 Pu Jfj]. However, very short lifetimes, in the range from 10 ps to 10 
ms, and experimental difficulties precluded, and in fact preclude to this very day, gainining 
sufficient knowledge on the nuclear structure in the second minimum. Only four values of the 
quadrupole moment in the second well were measured via nuclear methods of the total of 34 
shape-isomeric states. Three additional data for Am isotopes come from the optical isotope 
shift and hyperfine structure measurements of the difference between the mean square radii 
in the ground and isomeric states |7|. Moreover, difficulties in the experimental access to 
fission isomers make the results very dependable on rather extended argumentation chains 
involving assumptions instead of established experimental facts. Although this does not 
necessarily invalidates the claimed results, quite a number of legitimate reservations may be 
formulated which diminish their firmness, see js]. 

We determine properties of low-lying negative-parity shape vibrations in even-even ac- 
tinides using the collective model, i.e. the Schrodinger equation with deformation parameters 
as coordinates. As the ingredients we take the microscopic-macroscopic energy and cranking 
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mass parameters. In order to describe nonaxial modes we include nonaxial deformations of 
the nuclear shape corresponding to K — 1, 2, 3. The use of the microscopic-macroscopic 
energy is consistent with the working assumption that it correctly predicts the shape depen- 
dence of energy for "cold" nuclear configurations. On the other hand, using adiabatic mass 
parameters is only an approximation. It should work reasonably well for phonon energies 



he experimental energies in 
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smaller than 2 A, with A the neutron or proton pairing gap, 
the second well and many in the first well fulfil this condition 

One might expect that negative-parity shape vibrations are mostly octupole. However, 
there are at least two reasons to consider admixtures of higher order odd-multipolarity 
modes, i.e. A = 5 and 7: 1) The equilibrium shapes are spheroidal, hence the multipole 
components of different A are not the normal modes (not even being orthogonal) and may 
couple with each other, 2) Deformation parameters ftx^ of the microscopic-macroscopic model 
define the s.p. potential, not the density. Therefore, the octupole part of the latter may 
be produced by various odd-A deformations of the potential. As a consequence, one should 
account for a possible mode coupling when looking for low lying excitations of negative 
parity. We are not aware of any other study of shape oscillations explicitely including odd- 
rank multipoles of the order higher than three. 



One can observe that an alternative aj 
schematic Random Phase Approximation [9 



broach to the study of octupole vibrations, the 



while free from the adiabaticity assumption 



that we use for mass parameters, has its own problems, including the necessity of fixing the 
coupling constants. In order to study the coupling of various multipolarities, one would have 
many constants to fix. 

II. METHOD 

We use a sigle-particle (s.p.) Hamiltonian with the deformed Woods-Saxon potential 



defined in terms of the nuclear surface, according to the scheme exposed in [12J. However, the 
potential used in the present work admits more general nuclear shapes: The only restriction 
imposed on them is that they have one symmetry plane y-z. This leaves one conserved 
signature quantum number, s x = ±i, being the eigenvalue of the signature operator S x = 
VTZ^ 1 , with V the intrinsic parity and 1Z X the rotation by 7r about the intrinsic x axis. The 
degeneracy of pairs of states with s x = ±i (the Kramers degeneracy) reduces by half the 
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dimension of the s.p. Hamiltonian matrix. 

Nuclear shapes compatible with the assumed symmetry are defined by the following 
equation of the nuclear surface 

R(8,<p) = c({(3})Roa + Y,^oY X o(0,v) 

A>1 

+ E PxvY^e,?) 

\>l,/i>0,even 

+ E /V*W*.vO). (1) 

\>l,fi>0,odd 

where c({/3}) is the volume-fixing factor. The real-valued spherical harmonics, Y\ MC with 
even \x > and Y\ MS with odd /i > 0, are defined in terms of the usual ones as: Y^c = 
{Yxn + Yx-n)/\/2 and Y Ams = —iiY\^ + Y\_ M )/\/2. In other words, the dependence of the 
shape on the azimuthal angle ip enters through functions cos(fup) with /i even and sm(np) 
with fi odd. 

For the macroscopic part we used the Yukawa plus exponential model [lsj]. The param- 
eters of both the macroscopic part and the s.p. potential used in the present work, as well 
as the way the shell- and pairing corrections are calculated, are the same as in a number of 



previous studies, e.g. in [14]. 



Oscillations around I and II minima 



The second minima in actinides found in this and previous calculations, see e.g. [15|, Il6| . 
correspond to the axially- and reflection-symmetric shapes. The same holds for the g.s. 
minima, except for some light thorium, uranium and plutonium isotopes with N < 138 that 
have octupole equilibrium deformations and are not included in the present study. So we 
are left with nearly parabolic, not too shallow minima. Small oscillations around them with 
different intrinsic K numbers or parities are nearly independent. Indeed, the amplitude of 
vibrations cannot be too large, while exactly at the minimum the modes are uncoupled, 
except for the Coriolis coupling which may be then considered as a perturbation. 

The treatment of negative-parity shape oscillations in the collective model is based on 
the collective Hamiltonian 



V{Pk) (2) 
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diagonalized within the space of collective wave functions with the scalar product (ipi 



the mass tensor with indices corresponding to deformations /3\k, A = 3, 5, 7. In the present 
work we use cranking mass parameters. Oscillations for each K = 0,1,2,3 are considered 
separately. A similar approach was used in to study the K = octupole state in 222 Ra. 
In contrast to that work, where both quadrupole and octupole coordinate were used, we 
restrict collective variables to reflection-asymmetric deformations. 

Since our study is confined to nuclei with sufficiently deep, nearly parabolic minima, 
we can use the approximation of small oscillation amplitudes. This formalism follows from 
the one above if we replace V by the quadratic form (1/2) YUj CijPiPji approximately valid 
around the potential minimum, with the stiffness coefficients Cy, and fix the mass parame- 
ters at the values calculated at this minimum. Then the Hamiltonian becomes 



Within this 3D harmonic oscillator model, the study of the A = 3, 5, 7 coupling is straight- 
forward. The eigenmodes are given by a transformation of coordinates 



where the orthogonal matrices Si and S 2 diagonalize the mass tensor, SfBSi = B D , and 
the frequency matrix Q, SfYlS^ = Qd, with Q 2 = Bj^^ 2 SfCSiBj^ 1 ^ 2 . The square of the 
oscillation frequency Huk corresponding to the oscillation with a given K =0,1,2,3 is the 
smallest solution to the cubic equation det(Cjj — (huK) 2 Bij) = 0, with C^- and the 
stiffness and mass 3x3 matrices, i,j = (3\k and A = 3, 5, 7. 

The calculations of the coupled oscillations with the Hamiltonian fl2]) were made only for 
selected they require a time-consuming calculation of mass parameters, especially 

for if / 0. Some technical aspects of these calculations, as well as the evaluation of the 
cranking mass parameters for the rich deformation set Eq. (TjQ), are shortly described in 
the Appendix A. As discussed in Sect. Ill below, the results obtained with deformation- 
dependent mass parameters are not very different from those of the much simpler version 
with constant mass parameters ()3]). 



^2) = / d 3 /3kV det Bipiip2- Here V(/3k) is the microscopic-macroscopic energy and Bij((3k) is 




(3) 
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B. Electric dipole transitions 



Reduced probabilities of electromagnetic (EM) transitions between the rotational band 
built on the one-phonon state and the g.s. band can be calculated assuming the fixed 
structure of both the phonon and the collective rotor [181 ] . For an operator Ai of the 
multipolarity A one has 

B(X; K x = 0, I x -+ K 2 , I 2 ) = (2 - 5 K2 o)(hOXK 2 \ I 2 K 2 ) 2 \ (K 2 | M(X, K 2 ) \K X = Q) | 2 (5) 

with M(X, K 2 ) the intrinsic spherical component. For the negative-parity shape vibrations, 
the most prominent are dipole transitions, if not hindered by special reasons, as the octupole 
transitions are usually much weaker. It is especially true for shape isomers where only E0 



and El transitions are observed [22]. We have M(E1,0) = [3/(47r)] 1/2 D z and M(E1, 1) = 
-[3/(47r)] 1/2 (Ar + iD v )/y/2, where the dipole moment D = e(NY, p r p - Z£ n r n )/A Since 
the K — 1 phonon is the equal-weight combination of the equal-energy phonons in directions 
x and y, the intrinsic matrix element in Eq. (jSJ) equals — [3/(47r)] 1 / 2 D*, where D y = (Y\i s \ 
D y | g.s.), with symbolically denoted one-phonon state induced by the deformations Y\ ls , 
A = 3,5,7, Eq.©. 

Transition matrix elements in the intrinsic frame between the g.s. and the lowest excited 
state of negative parity | it—) could be calculated by integrating the transition density 
PtriPk) = {Vdet Bipl_ip gs }((3k) with the proper operator, represented in the collective space. 
Taking the dipole moment as an example, we have 

D l = <tt- | D I 0} = / d 3 l3 k p tr (l3 k )D(/3 k ), (6) 

J — oo 

where D{j3k) is the expectation value of D in the mean field state with the deformations (3 k 



A considerably simpler approximation consists in calculating the diagonal matrix element 
of the transition operator in the mean-field state with the deformations (5^ fixed as the most 
probable by the above transition density: 



/oo 
d 3 ^ jPtr (^ k . (7) 
-oo 

The dipole and octupole moment operators are to the leading order linear in the odd- 
multipole deformation parameters (3\k, so indeed, the integration in Eq. (jSJ) replaces (3\k in 
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D({3zk, fibK, 07K, •••) by {3\ r K . A nonlinearity of D as a function of deformations may introduce 
an error in this approximation. Assuming that the lowest negative-parity mode corresponds 
to the normal coordinate £1, the values of /3|. r for harmonic vibrations are actually given by 
Eq.® with £f = l/y/2, 6 = 6 = 0. 

This may be contrasted with the strong coupling limit with two octupole minima, at 
±Pl 9 K , in which the transition matrix element D f is calculated as the expectation value at 
this deformation of equilibrium. One can notice that in this case, [3 e x q K is equal to ftxKi 
the expectation value of (3\k in the first excited state of negative parity, nearly degenerate 



with the g.s., (3^ = 2 / °° ydet B \ i/j^- 1 2 fi\K<$fik- For the harmonic lowest-lying phonon 
one has the relation ^ K = 0.63/3^ which follows from Eq. fll]) and the relation for the 
one- dimensional harmonic oscillator: £\ r = (vr/8) 1//2 ^ _ ~ 0.63^", as = 2/y/n. 

Expectation value of the electric dipole moment in the state with deformations 0\k, 



K =0,1, is calculated as a sum of the macroscopic and shell-corection parts, see e.g. [19l. |20] . 
The macroscopic part, derived within the Droplet Model in 21], has to be calculated as in 
22|], i.e. without assuming small Px^. For the completeness of the presentation, the relevant 
formulas and parameters are collected in the Appendix B. 



C. Estimate of the Coriolis coupling effect 

Within the rigid rotor-vibration coupling model, there are two modifications of the energy 
of stationary states with the total angular momentum I relative to the energy of one-phonon 
states: (i) a shift by a(I(I + 1) — K 2 ), with a = h 2 /(2<J), J the moment of inertia, and (ii) 
the effect of the Coriolis coupling 

(JM; K + l\ Hear I IM- K) = -ay/l + 6 Koy /(I -K)(I + K + 1)(K + 1 | J + | K) , (8) 

connecting states differing by AK =1. In the spherical limit, the matrix elements of J + 
between substates K and K + 1 of a collective vibration of the multipolarity A are equal to 
^(A — K)(X + K + 1). The actual matrix elements for octupole phonons in deformed nuclei 
were found close to the spherical limit in Here we also include multipoles A = 5 and 7, 
so we estimate the effect of the Coriolis coupling using the spherical value with A consistent 
with the composition of the lowest-lying phonon, see sect. III. 

For a test case of / = 1 states (the mixing of K = and K — 1) in 240 Pu, in the first 
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minimmum, we take A = 3 and with a = 7.156 keV 23| find the coupling if<7 or =49.6 keV; 
in the second well, with a = 3.33 keV jl], we obtain 7^^=23.1 keV for A = 3 and 36.5 keV 
for A = 5. As the energy shifts of I = 1 states due to the coupling (|S]) are smaller then 
the coupling itself, especially when the difference in Coriolis-unperturbed energies is much 
larger than Hc or , which is often the case, we do not include them in presented results. 



III. RESULTS AND DISCUSSION 

We begin with the description of energy minima obtained by a minimization in a multidi- 
mensional deformation space. A typical plot of energy vs. deformation fco after minimiza- 
tion with respect to other axial deformations is shown in Fig. [I] for 240 Pu. It is normalized 
by setting the macroscopic energy to zero at the spherical shape. The experimental ex- 
citation energy of the superdeformed state, determined by the statistical analysis, equals 
j^exp _ 2 MeV whereas our calculation gives E% = 2.0 MeV. The calculated first 
and second barrier, E 1 ^ = 6.4 MeV and E 1 ^ = 5.0 MeV, are in a satisfactory agreement 
with the experimental values E e £ v = 5.8 MeV and E e ^ p = 5.45. As mentioned before, 
mass-asymmetric and nonaxial deformations do not change the first and second minima. 

One can also see in Fig. 1 the calculated third minimum. A detailed discussion of 
the third minima goes beyond the scope of the present paper. Here, we can mention that 
nonaxial deformations, not included in Fig. 1, modify the energy surface around them. 
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FIG. 1: Calculated fission barrier for 240 Pu as a function of $ 
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Energies of shape vibrations depend on the values of the mass and stiffness parameters 
in the vicinity of the equilibrium deformation. The stiffness coefficients were determined by 
fitting energy by a quadratic form in deformations. A quality of this approximation in the 
first and second well is illustrated in FigfS]for the octupole deformation /3 30 . It can be seen 
that the assumption of harmonicity is satisfied in the reegion of excitation energies up to 
2-2.5 MeV. The root-mean-square deviation of the fit for 234 U equals 1.3 keV and 0.45 keV at 
the first and second minimum, respectively. Errors of the fit, for all nuclei and deformations, 
fall in the range 1-20 keV, with a typical number of points in each fit equal to 30. 
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FIG. 2: The extraction of the stiffness C3030 at the I-st and Il-nd minimum for U. 

A collection of all calculated K = stiffness coefficients is shown in Fig. [3] for the first 
minima and in Fig. @]for the second minima. There are diagonal C3030, C50505 C7070 and off- 
diagonal C3050, C307O) C5070 stiffness parameters. It is natural when the latter are smaller 
than the former. Very telling is the relative smallness of the C3030 coefficients at the first 
minima, where they are similar in magnitude to the off-diagonal coefficients, and their much 
larger values at the second minima. 

The plots of mass parameters -B3030 and -B5050 around the reflection- symmetric I and II 
minimum are shown in Fig. The 240 Pu nucleus has been chosen as an example. One can 
see that the assumption of constant mass parameters, fixed at their values at the minima, 
is not drastically wrong. There are two reasons why a deformation dependence of mass 
parameters should not be crucial: 1) The relevant deformation range is limited, as the peak 
of the probability density of the one-phonon state corresponds to (3^ < 0.1; 2) Phonon 
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FIG. 3: Stiffness coefficients at the I-st minimum. 
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FIG. 4: Stiffness coefficients at the Il-nd minimum. 

energy depends on the mass parameter as (C/B) 1 ^ 2 , which is not a very strong dependence. 
One can observe a different behavior of -B3030 around the I-st and Il-nd minima. While at 
the first minimum this parameter takes the maximal (or close to the maximal) value, it 
reaches the minimal value at the second minimum. Both -B3030 and -B5050 are skewed with 
respect to the /3 30 and /3 50 axes. 

Results of the simpler, harmonic model (j3]) that already contain the effect of the coupling 
of various multipolarities are presented below. Their modification due to the variability of 
the mass parameters was worked out for selected cases and is discussed afterwards. 
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FIG. 5: Cranking mass parameters -B3030 , -B5050 around the first (top panels) and second (bottom 
panels) minimum for 240 Pu in the ( /33o,/35o) plane. 

A. Results with constant mass parameters 

Calculated energies of the lowest negative-parity shape oscillations in the first minimum 



for various even-even actinides are shown in Fig. [6] together with experimental results 24]. 
Only vibrational cases are considered which means that we excluded nuclei with the octupole 
equilibrum deformation or very close to the reflection-symmetry breaking. As the calculated 
energies of the K = 3 phonons are greater than 1.5 MeV in Th, U, Pu and Cm isotopes and 
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greater than 1.3 MeV in Cf isotopes, and, moreover, an experimental information on such 
octupole states is scarce and uncertain, we do not show these results here. It may be seen 
that for the K = and K = 1 modes an overall agreement to within ~100 keV (except for 
230 U) is obtained between the data and theoretical evaluation. It is worth mentioning that 
our calculations do not involve any adjustment to the data they are compared with. For the 
K = 2 mode the experimental energies are considerably overestimated. On the other hand, 
one has to notice that the experimental assignments made for K = 1 and K = 2 octupole 
band-heads are not many and some of them seem quite uncertain. One can observe in Fig. 
E]a clear change in energies between iV = 146 and N = 148(150): they increase for K w = 0~ 
and 1~ in U, Pu, Cm and Cf, and decrease for K 71 = 2~ in U, Pu and Cm. 

From the difference between the energies and the quantities {C^kzk / B^kzk) 1 ^ 2 one can 
judge the importance of the off-diagonal components of the stiffness and mass tensors. The 
resulting mixed multipolarity character of the lowest phonon may be also seen from the 
deformation parameters (3 1 { K that follow from Eq.(j4]). For the K = mode, they show a 
small admixture of ft^o and/or ft-jo deformation to /?3o: A = 5 mixing is important for the 
lightest/heaviest isotopes, A = 7 for the middle. For example, we have (/3|q, j5f Q , /3|q) = 
(0.049, 0.013, -0.015) in 232 U, (0.058, 0.005, -0.017) in 238 U and similar values in 240 Pu and 
246 Cm, while the fo -kv coupling is implied by the values (0.052,-0.017,0) in 250 Cf. For 
the K — 1 mode, the admixture of ft^i is, except for Th isotopes, larger than that of /5 71 : 
in 238 U we have /3j^ = (0.050, —0.016, 0.007); similar values occur for 240 Pu and many other 
nuclei. As shown in Fig. [6j the role of the coupling to higher multipolarities decreases with 
the increasing i^-number in the first potential well. In particular, for K w = 2~ vibrations, 
one could use the simple formula hu>2 = f° r the excitation energy. Appreciable 

admixtures, mainly of A = 5, are present for the lighter isotopes. 

In Fig. [7J calculated low-lying negative-parity excitation energies at the second minimum 
are shown as a function of the neutron number for different isotopes. The plotted experi- 
mental data were taken from [l|. The K = 3 energies are not presented as they lie above 
2.5 MeV, even higher than in the first minima. We do show calculated vibration energies 
in the second well in Cf isotopes, although no shape isomer was established there up to 
now. In our calculations, the barrier protecting the second minimum in 246 Cf against fission 
is by 1 MeV smaller and thinner than in its isotone 244 Cm, which supports the short-lived 
(< 5 ps) shape isomer and its longer lived (< 100 ns), probably spin-isomeric, 1.3 MeV 
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excitation 



251 ] . However, it is not completely excluded that, with a better technique, some 



shape isomer in Cf could be established in the future, at least in an odd isotope. 

Rather extensive experimental results were collected on the excitations in the II minimum 
of 240 Pu j^fl [27]]. These include negative-parity rotational bands as well as the candidate for 
the /3- vibration band. The lowest rotational bands interpreted as K 71 = 0~, 1~~ and 2~ have 
bandheads at the excitation energy 555 (I 71 = 836 (J 71 " = 1~) and 806 keV (I n = 2~), 
respectively. Only two first states, I n = 1~ and 3~, of the supposed octupole K n = 0~ band 
were observed. There are three more bands identified as K n = l - : at 936, 1246 and 1344 
keV; the supposed /3 band starts with the second I n = + state at 770 keV. 

As follows form FigJTJ within the presented model we cannot reproduce the K = 0~ 
excitations at the second minimum. In 240 Pu, we obtain energy nearly three times too 
large. Situation is even worse for 236 U, in which the back-decay to the first minimum was 
experimentally detected and the excitation energy of mere 290 keV was attributed to the 
K n = 0~ state in the second well. One reason for our results can be seen in Figs. [3] and 
IH by comparing the stiffness C3030 against the octupole deformation /3 30 at the first and 
second minima. For all studied nuclei, the stiffness in the second minimum is roughly two 
times larger than in the first one and mostly larger than the stiffnes C5050 against P^. The 
second reason for the larger K = excitation energies in the second well comes from the 
mass parameters: -B3030 are two to three times smaller at the second minimum than at the 
g.s., and even little smaller than -B5050? opposite to the situation at the first minimum, see 
Fig. El 

As a result, the lowest K n = 0~ excitation in the second well is a mixture of the /^o and 
f3$o modes, with the latter dominating for nearly all isotopes. For example, the values of 
are equal to (0.033,0.040,0) in 230 Th, (0.032,0.041,-0.005) in 236 U, (0.030,0.040,-0.007) 
in 240 Pu, (0.024,0.041,0) in 248 Cm and (0.020,0.040,0.003) in 250 Cf. When the coupling 
between different multipolarities is neglected, the /3 50 mode is the lowest in the second 
minimum. 

There is a large energy lowering below the uncoupled value (C303o/-B303o)^ 2 due to this 
coupling, see Fig. [71 This shows that it is impossible to describe low-lying K n = 0~ states 
in the second well as a pure octupole vibration. However, the coupling with the Y 50 mode 
is not sufficient in view of the experimental results. 

The couplings for the K 11 = 1~ and 2 _ modes reduce phonon energies by ~ 300 and 
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~ 400 keV, respectively, below the values of (C 3 K3K I Bzkzk) 1 ^ 2 > Fig. [71 This is much less 
than for the 0~ mode. The values P f { K show a conspicuous mixing of A = 3 and 5, with much 
smaller admixtures of A = 7. For K 71 = 1~ phonons, the octupole mode always dominates, 
/^3i > Aji- m m °st cases this is also true for the K 71 = 2~ mode, but for some nuclei 
reaches the size of Both phonon energies in 240 Pu are overestimated, by 400 (l - ) and 
170 keV (2 _ ); the supposed 2~ band head at 830 keV in 236 U is also calculated too high 
by 430 keV. The difference in the calculated K 71 = 2~ phonon energies in 236 U and 240 Pu 
follows mainly from the larger mass parameter -B3232 in the latter. It is remarkable that in 
the second minimum the mass parameters -B3131 > -B3232 are ~2 or more times larger than 
S3030, while the latter are the largest in the first minimum. 



B. Results with shape-dependent mass parameters 

Oscillation energies Huk in the first and second minima calculated with shape-dependent 
mass parameters for selected even-even actinides are compared in Table |T] to the constant- 
mass variant of the calculations and the experimental data, where available. As may be seen, 
the model including the shape dependence of the mass parameters (T5]) spoils or improves 
the agreement with the data, depending on a particular case. In the first well, there is a 
substantial decrease in some K n = 0~ energies away from the measured values, especially for 
240 Pu. Some other energies become closer to the data, e.g. in 236 U and 246 Cm. The K 71 = 1" 
energy in 230 Th is spoiled, while for the = T~ mode there is an improvement for 240 Pu 
and worsening for 250 Cf, already 400 keV off the data in the harmonic approximation. These 
changes follow from the differences between the mass parameters at the minimum and their 
averages over a region around it. 

In the second minimum, K n = 0~ energies from the full Hamiltonian fl2]) are by 200-300 
keV lower than for its simplified version (j3j), except for 234 U, for which Uujq stays the same. 
This is an improvement, but not a sufficient one, as now the = 0~ energy for 240 Pu is 
roughly 2.4 times too large instead of three. The = 1" and 2~ energies in 236 U and 240 Pu 
do not change much, except for the K n = 2~ state in 240 Pu, that increases further 200 keV 
away from the experimental value. 

The structure of vibration phonons changes mainly according to the energy change. The 
values of /3^ K increase appreciably with respect to the constant mass parameters version only 
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TABLE I: Energies of negative-parity shape vibrations (in keV) for K = 0,1,2 around I and II 
minimum obtained from the diagonalization of the full Hamiltonian ([2]) - H co n, using constant 
mass parameters in the minimum (E2) - 

Experimental data are taken from [2j] for the I-st minima and from [l| for the Il-nd minima. 



Bconst and calculated by a simple formula: Tiujk 
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K = 
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92 
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527 


756 


367 


1138 


1309 


2482 




92 


232 


698 


546 


773 


563 


1132 


1331 


2430 




92 


234 


820 


689 


906 


786 


1400 


1390 


2672 




92 


236 


694 


644 


881 


688 


1199 


1430 


3002 


290 


92 


238 


451 


561 


795 


680 


1236 


1492 


3104 




94 


238 


465 


529 


814 


605 


1197 


1505 


2999 




94 


240 


338 


509 


792 


597 


1272 


1565 


3254 


555 


94 


244 


1062 


1083 


1129 


950(?) 


1353 


1654 


3339 




96 


246 


1092 


1138 


1217 


1079 


1380 


1681 


3440 






K = 1 


K = 1 


90 


230 


770 


915 


1033 


952 




962 


1373 




92 


236 




1036 


1175 


967(?) 


1354 


1461 


1765 




94 


240 


858 


859 


1015 


938(?) 


1266 


1289 


1568 


836 




K = 2 


K = 2 


90 


230 


1264 


1253 


1376 


1079(?) 




1367 


1760 




92 


234 


1377 


1348 


1488 


989 




1287 


1704 




92 


236 




1385 


1490 


1110(?) 


1345 


1226 


1639 


830(?) 


94 


240 


1295 


1469 


1497 


1241(?) 


1166 


972 


1355 


806 


98 


250 


1303 


1274 


1277 


872 




1061 


1738 





if there is an appreciable decrease in the phonon energy. Their ratios are similar as those for 
the harmonic Hamiltonian. It should be noticed that for anharmonic vibrations the simple 
relation between f3^ K and the multipolarity composition of the vibration coordinate, based 
on Eq.Q, is lost. 

One can say that the results obtained by the full diagonalization of the Hamiltonian ([2]) 
do not introduce drastic changes in the results obtained with ([3]). 



C. Electromagnetic dipole transitions 

Transition dipole moments governing El transitions between one-phonon K = and 
K = 1 bands and the g.s. band have been calculated using the values f3^ K from the model 
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(J3]) with constatnt mass parameters. They are given in Table [TT] for selected nuclei. At 
the first minimum, D* values show a considerable variation from a nucleus to nucleus. The 
K — and 1 components are correlated, e.g. both are large for 230 U; this large K = value 
is roughly consistent with the measurements [28]. For K — 1, B(E1) ~ 2(D t ) 2 , Eq. (J5J), 
so for an easy comparison of -B(ET)s for if = and K = 1, one should multiply £)* for 
if = 1 from Table El by v^- For some nuclei, the opposite signs of the microscopic and 
macroscopic parts lead to a nearly complete cancellation of D l . This results in a variation 
of the deexcitation pattern along the isotopic chain, as may be seen for 23 °- 238 u. 

More specifically, the negative-parity rotational band decays according to the ratios of 
the rates of E2 transitions within, and El transitions out of the band. From the calculated 
transition dipole moments D t in [e fm] we can estimate the ratio 

T(E1),^ _ L303 |||l) ,8, (2/- l)(/- l) (VV _ 



T(E2) M - 2 E'-(E2) V57 (I-lf-IC- \Q„ 

where Qo is the intrinsic quadrupole moment of the negative-parity band in units of [10 b], 
and energies E 7 of gamma transitions are given in MeV. Using values of Qo measured for 
the g.s. bands, we obtain the following T(E1)/T(E2) ratios for transitions from the I 71 = 7~ 
state of the K w = 0" band: around 300 in 232 > 234 U, Lg in 236tj and arounc i 10 o in 238 U. 
Although not a perfect match with tha data, this roughly correlates with no intraband E2s 
seen along the negative-parity bands in 232 ' 234 U ; fa e complete regular E2 cascade and no Els 
above the 3~ state in 236 U, and the E2 cascade ending at the 7~ state, with more Els in 238 U 



29 



31] . On the other hand, the predicted nearly perfect cancellation of D l for K = in 



240 Pu is not supported by the value ?»0.12 efm measured at spin 11 and the smaller values 
for the neighbouring Pu isotopes [32]. The relatively small dipole moment measured for 



230 Th [33] does not contradict our K = 1 value. 

In the isomeric second minimum the situation is entirely different. The macroscopic part, 
Dmac ~ faofizK in the leading order, is dominant there as its magnitude becomes larger than 
that of the shell-correction part due to the large equilibrium value of /?2o- F° r K = bands, 
the microscopic part adds to the macroscopic part. The large values of the dipole moment 
and small energies of rotational transitions near the bandheads preclude observation of the 
intraband E2s in negative-parity vibrational bands in the isomeric minima in actinides. This 
agrees with the lack of experimentally observed E2 transitions in the isomeric minimum in 

24 °pu Q,y,0. 
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TABLE II: Macroscopic D 

mac-, microscopic D m ic and total transition dipole moment Eft for K v — 



and 1 phonons in the first and second minimum. 
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92 
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0.19 


0.04 


0.23 


0.29 


0.11 


0.40 


92 


232 


0.17 


0.00 


0.17 


0.34 


0.08 


0.42 


92 


234 


0.17 


-0.07 


0.10 


0.34 


0.04 


0.38 


92 


236 


0.17 


-0.16 


0.01 


0.34 


0.02 


0.35 


92 


238 


0.17 


-0.24 


-0.07 


0.33 


0.02 


0.35 


94 


238 


0.19 


-0.13 


0.06 


0.33 


0.09 


0.42 


94 


240 


0.18 


-0.18 


0.00 


0.32 


0.08 


0.40 


94 


244 


0.15 


-0.18 


-0.03 


0.29 


0.03 


0.32 


96 


246 


0.15 


-0.14 


0.01 


0.29 


0.12 


0.41 




K = 1 


K = 1 


90 


230 


-0.042 


0.031 


-0.01 


-0.157 


-0.077 


-0.23 


92 


230 


-0.049 


-0.079 


-0.13 


-0.173 


-0.055 


-0.23 


92 


232 


-0.047 


-0.041 


-0.09 


-0.179 


-0.035 


-0.21 


92 


234 


-0.049 


0.008 


-0.04 


-0.174 


-0.017 


-0.19 


92 


236 


-0.055 


0.035 


-0.02 


-0.170 


0.00 


-0.17 


92 


238 


-0.063 


0.054 


-0.01 


-0.174 


0.00 


-0.17 


94 


238 


-0.059 


0.024 


-0.036 


-0.182 


0.023 


-0.16 


94 


240 


-0.065 


0.045 


-0.02 


-0.183 


0.022 


-0.16 



IV. CONCLUSIONS 

Our study of low-lying negative-parity oscillations in even-even actinides leads to the 
following conclusions: 

• Considering that we have no adjustable parameters, the data on negative-parity exci- 
tations in the first minima are reproduced quite well for the K w = and 1~ phonons; 
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energies of the 2 phonons are significantly overestimated. 

The model predicts = 0~ energies in the second well in 240 Pu and 236 U that are 
three and more times larger than the claimed experimental values. For the K =1,2 



ll) for 



phonons, the agreement with the data is better, similar to that obtained in 
240 Pu, but still the calculated energies are too high. 

• Low-lying negative-parity oscillations show small admixtures of the multipolar it ies A = 
5,7 to the octupole mode at the first minima, and equal or dominant contribution of the 
A = 5 multipole at the second minimum. Hence the coupling of various multipolarities 
is important in the description of the "octupole" vibrations, especially in the second 
minimum. 

• Taking into account the multipoles A =3,5,7 in the phonon structure and the exact 
macroscopic contribution to the dipole moment, we predict large transition dipole 
moments from the "octupole" band to the g.s. band in lighter actinides at the first 
minima, and for all investigated nuclei at the isomeric minima. 

We do not see any way to reconcile our model with the reported data on the energies 
of the K n = 0~ mode in the second well. Both the calculated sizable stiffness Fig. H] and 
small mass parameters Fig. [5] suggest that either our model is completely unreliable there, 
or, perhaps, the experimental K n = 0~ assignments in the shape isomeric minima in 236 U 
and 240 Pu are not related to the collective shape vibrations. 

One could think of possible improvements of the model. A natural step would be to 
include the quadrupole-octupole coupling. However, as long as the second minima are 
reflection-symmetric, this would be a second-order effect, while the energy surfaces and 
mass parameters do not hint to its unusual enhancement. One could also consider a fine- 
luning of the pairing strength in the second well, or including the quadrupole pairing as in 



111 ] . Still, the required reduction of the 0~ phonon energies is so large, that the pairing alone 
hardly can be a cause. On the other hand, the discrepancies observed for 1~ and 2~ phonons 
in the Il-nd well probably could be reduced by a fine tuning of the model parameters. 
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Appendix A: Some aspects of the calculations 



The matrix diagonalization of the Hamiltonian fl2]) is performed in the basis of the three- 
dimensional harmonic oscillator. The matrix elements of the Hamiltonian were calculated 
with the 20 point Gauss-Hermitte quadratures. The calculations were reduced by setting 
the mass parameters to their values at the minimum outside the mesh of 14 x 14 x 14 points 
for the easier K = case, and outside the mesh of 8 x 8 x 6 points for K = 1 and 2. We have 
checked that the related distortion of the Hamiltonian is unimportant for the lowest phonon 
states, as even the smaller mesh covers their peak region. As a check of the program, we 
reproduced the harmonic case of constant mass parameters. We also checked the case of the 
position-dependent mass parameter B ~ (3 2 with the quartic potential ~ /3 4 , which gives 
the harmonic spectrum. 

The cranking mass parameters were calculated by replacing derivatives with respect to 
deformations by finite differences. This has some advantages for many deformation pa- 
rameters, but introduces an error which may reach 2-3% for diagonal components. The 
new approach was tested with the old code for mass parameters for axial and some non- 
axial deformations. The main tool of these calculations, the code that diagonalizes the s.p. 
Hamiltonian, was checked independently with the older, less general versions. 

An independent test of this code is provided by the fact, that different deformation 
sets may correspond to the same shape, so they should produce the same s.p. spectrum. 
In particular, the deformations Y\Ks with K = 1,3, X = 3,5,7 of an axially symmetric 
equilibrium shape give the same s.p. spectrum as the deformations Y\k c - The latter are 
accomodated by the parametrization Eq. ([T]), by switching the choice of the symmetry axis 
from z to y. This follows from the relations between spherical harmonics defined with respect 
to the reference axes z, x, y (denoted as Y y ) and the standard ones: 
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5V7 15V7 3v/33 



v . 3v/2T 19 3v^43 

73c " ~^2 70 ~ ^2 72c + "32" F74c + 76C ' 



The equivalent deformation set for axially symmetric first and second minima follows from 
the expressions for spherical harmonics Y^ : 

^40 = g^40 + -^-Y42c + — g— ^44c, (A2) 

y _ 5 v/105 3v/7 v/23T 



35 3^35~ 3a/77 429 Sv^, 

128 80 + 3272 82c + ~MT Y84c + 3^2 86c + ~T28~ 



Appendix B: Dipole moments within the microscopic-macroscopic method 

The macroscopic part of the expectation value of the electric diple moment is calculated 



as a sum of the redistribution and the neutron skin effects 



2l| 



D — D re d + Dsfcin; (Bl) 

where for Z protons and N neutrons, A = Z + N, I = (N — Z)/A, one obtains from the 
Droplet Model 

AZe 2 (\ 6LI\ 



D 



red 



J KJ 



j (< v > v < £ >v - < < >v) , (B2) 



2NZ 

Vskin = - ~5)Rq (< £ > v - < £ >s) (B3) 

9 ZA 2 ' ?J e 2 

+ q Bs (< v >s< ^ >s - < < >s) • 



The Droplet Model parameters are: the volume symmetry-energy coefficient J, the nuclear 
incompressibility K, the effective neutron skin stiffness Q, the density symmetry coefficient 
L, the nuclear radius Ro = r^A 1 ^ and the equilibrium value of the average relative neutron 

excess, 

_ i + *?ZA- 2 I* 

*= rffe ■ (B4) 
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Nuclear shape enters through the constant Bs, being the ratio of the area of a deformed 
surface to the surface area of the sphere of the same volume, and the averages, over both 
the nuclear volume and surface, of the scaled radius vector £ = t/Rq and the Coulomb 
potential v in units of Ze/Ro, with < f >v meaning j v f/V, and < f >s meaning J s f /S. 
The calculations were performed with the following values of the parameters: r = 1.16 fm, 
J =32.5 MeV, K =240 MeV, Q =50 MeV and L =100 MeV, as in [22J. 

The microscopic part of the dipole moment is calculated as D micr = C(< D > —D) 

nn 

[191 . 1201 ] , where the first part is the expectation value of the dipole operator on the deformed 
state and D is the analogous expression in which the actual pairing occupation numbers 
2vf were replaced by the quantities smoothed according to the Strutinsky prescription. The 
factor ( takes care of the reduction of the effective charge due to the particle-vibration 



coupling to the E 



calculations 



20|,|22|. 



giant resonance. We have used ( = 0.33 as it was done in previous 
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FIG. 6: Energies of negative-parity shape oscillations in the g.s. minimum for K =0,1,2. Results 
including all (ALL) multipolarities A = 3,5,7 are denoted by squares; triangles follow the formula: 
Tiuj k = ^ bIkIk ' Ex P erimental data (EXP) 



uncertainty in the experimental assignment. 



24j ] are marked with circles; parenthesis signals some 
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FIG. 7: The same as in Fig[6l but in the second minimum. Experimental data (EXP) are taken 
from 
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